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We consider minisuperspace models constituted of Bianchi I geometries with a free massless 
scalar field. The classical solutions are always singular (with the trivial exception of flat space- 
time), and always anisotropic once they begin anisotropic. When quantizing the system, we obtain 
the Wheeler- De Witt equation as a four-dimensional massless Klein-Gordon equation. We show that 
there are plenty of quantum states whose corresponding bohmian trajectories may be non-singular 
and/or presenting large isotropic phases, even if they begin anisotropic, due to quantum gravitational 
effects. As a specific example, we exhibit field plots of bohmian trajectories for the case of gaussian 
superpositions of plane wave solutions of the Wheeler-DeWitt equation which have those properties. 
These conclusions are valid even in the absence of the scalar field. 
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I. INTRODUCTION 

Two of the main questions one might ask in cosmology are: 

1) Is the Universe eternal or it had a beginning, and in the last case, was this beginning given by an initial 
singularity? 

2) Why the Universe we live in is remarkably homogeneous and isotropic, with very small deviations from this 
highly symmetric state? 

The answer given by classical General Relativity (GR) to the first question, indicated by the singularity theorems 
0, asserts that probably the Universe had a singular beginning. As singularities are out of the scope of any physical 
theory, this answer invalidates any description of the very beginning of the Universe in physical terms. One might 
think that GR and/or any other matter field theory must be changed under the extreme situations of very high energy 
density and curvature near the singularity, rendering the physical assumptions of the singularity theorems invalid near 
this point. One good point of view (which is not the only one) is to think that quantum gravitational effects become 
important under these extreme conditions. We should then construct a quantum theory of gravitation and apply it 
to cosmology. For instance, in the euclidean quantum gravity approach Q to quantum cosmology, a second answer to 
the first question comes out: the Universe may have had a non-singular birth given by the beginning of time through 
a change of signature. 

In the same way, the naive answer of GR and the standard cosmological scenario to the second question is not at 
all satisfactory: the reason for the Universe be highly homogeneous and isotropic is a matter of initial conditions. 
However, solutions of Einstein's equations with this symmetry are of measure zero; so, why the Universe is not 
inhomogeneous and/or anisotropic? Inflation is an idea that tries to explain this fact. Nevertheless, in order for 
inflation to happen, some special initial conditions are still necessary, although much more less stringent than in the 
case without inflation. Once again, quantum cosmology can help in this matter by providing the physical reasons for 
having the initial conditions for inflation. 

Recently, many papers have been devoted to the application of the Bohm-de Broglic interpretation of quantum 
mechanics [|]|]] to quantum cosmology ■ One of the aims of these papers was to try to give an answer to the first 
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question within a new perspective in quantum cosmology. In the framework of the Bohm-de Broglie interpretation of 
quantum mechanics applied to minisuperspace models, it was shown in References [[ll] [l5| that quantum gravitational 
effects may indeed become important under extreme situations of very high energy density and curvature, and avoid 
the initial singularity of the classical models. The solutions that come out are then non-singular (eternal) cosmological 
models, but with a very hot phase, tending to the usual classical solutions when the energy scales become smaller 
than the Planck scale. Of course this is a much better answer than the one given by classical GR. 

The aim of this paper, as a natural sequence of the references cited above, is to address the second question 
from the point of view of the Bohm-de Broglie interpretation of quantum cosmology. We want to investigate if 
quantum mechanical effects can isotropize an anisotropic cosmological model yielding an alternative mechanism for 
the isotropization of the Universe. Our strategy is to take classical anisotropic models which never isotropizes, quantize 
them, and examine if quantum effects can at least create isotropic phases. As we will see, quantum effects can not 
only create isotropic phases but they can also isotropize the models forever, for a large variety of quantum states. 
Other authors have also studied this problem, mainly taking Bianchi IX models, adopting other interpretations of 
quantum cosmology, and they arrive at similar results [l7|-pl|]. 

The mathematical models we take are Bianchi I models with and without matter. For the case without matter, 
one of the anisotropic degrees of freedom is suppressed and we end with a two-dimensional minisuperspace model. 
The classical solutions of this case are all singular and always anisotropic, with the trivial exception of flat space- 
time. When quantizing this system we find that the solutions of the Wheeler-DeWitt can be written in terms of 
superpositions of plane wave solutions. Taking gaussian superpositions we show that there are bohmian trajectories 
representing universes expanding from an initial singularity and others representing bouncing non-singular universes, 
both with isotropic phases during the course of their evolutions. There are also bohmian non-singular periodic 
solutions with periodic isotropic phases and non-singular bouncing solutions which are never isotropic. For the case 
with matter, we take a minimally coupled scalar field with arbitrary coupling constant uj. The minisuperspace is 
four-dimensional. The classical solutions are again singular, and once they begin anisotropic they continue to be 
anisotropic forever. When quantizing the system, we arrive at a Wheeler-DeWitt equation which corresponds to a 
massless Klein-Gordon equation in four dimensions. We show that there are plenty of spherical wave solutions whose 
bohmian trajectories represent expanding universes that isotropizes permanently after some period, many of them 
which are also non-singular, gaussian superpositions also present this feature. 

The paper is organized as follows: in the next section we present the minisuperspace models which we will investigate 
and their classical solutions. In section III we quantize these models, and obtain their corresponding Wheeler-DeWitt 
equations and general quantum solutions. In section IV we present the results concerning the quantum isotropization 
of those solutions through the investigation of the bohmian trajectories we obtain from them. We end up in section 
V with conclusions and discussions. 



II. THE CLASSICAL MINISUPERSPACE MODEL 



Let us take the lagrangian 



R — wc 



(2.1) 



For w = — 1 we have effective string theory without the Kalb-Ramond field. For w = —3/2 we have a conformally 
coupled scalar field. Performing the conformal transformation g^ u — e^g^ we obtain 



L= V=9 



R — Cm, 



(2.2) 



where the bars have been omitted, and C w = (w + §). The cases of interest correspond to C w > 0. 

The gravitational part of the minisuperspace model we study in this paper is given by the homogeneous and 
anisotropic Bianchi I line clement 



ds 2 = -~N 2 {t)dt 2 + exp[2/3 (U + 2/3+ (i) + 2>/3 / 9_(*)] dx2 
exp[2#)(t) + 20+(t) - 2>/3/3-(<)] dy 2 + exp[2/3 (i) 



4/3+ (t)} dz 2 



(2.3) 



This line element will be isotropic if and only if /3+(i) and (3-{t) are constants (!]]. Inserting Equation (2.3) into the 
action S = J Ld A x, supposing that the scalar field <j> depends only on time, discarding surface terms, and performing 
a Legendre transformation, we obtain the following minisuperspace classical hamiltonian 
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H =2i£m { - pl+p2++p - +plh {2A) 

where (po,P+,P-,P0) & r e canonically conjugate to (/3o, /3+, f3— , <ft), respectively, and we made the trivial redefinition 

We can write this hamiltonian in a compact form by defining = (/?□, /?+, 4>) an d their canonical momenta 
Pn = (Po,P+,P-,P<f>), obtaining 

where rf ,v is the Minkowski metric with signature ( — V H — V). The equations of motion are the constraint equation 
obtained by varying the hamiltonian with respect to the lapse function iV 



and the Hamilton's equations 



U = if v p»p v = , (2.6) 

*'=S=Tdk>* /> *' (2 - 7) 
= -w = ■ (2 - 8) 

The solution to these equations in the gauge N = 12exp(3y ) is 

2/" = rf v Vvt + C , (2.9) 

where the momenta p v are constants due to the equations of motion and the are integration constants. We can 
see that the only way to obtain isotropy in these solutions is by making p\ = p + = and pi = p- = 0, which 
yields solutions that are always isotropic, the usual Friedmann-Robertson- Walker (FRW) solutions with a scalar field. 
Hence, there is no anisotropic solution in this model which can classically becomes isotropic during the course of its 
evolution. Once anisotropic, always anisotropic. If we suppress the <fi degree of freedom, the unique isotropic solution 



is flat space-time because in this case the constraint (2.6) enforces po = 



To discuss the appearance of singularities, we need the Weyl square tensor W 2 = W alillu W a p^ v . It reads 

W 2 = ^e- 12 ^(2p oP l ~ 6p p 2 _ P+ +p 4 _ + 2p\p 2 _ +pX+p 2 P 2 + + p 2 p 2 _) . (2.10) 



Hence, the Weyl square tensor is proportional to exp (— 12/3 ) because the p's are constants (see Equations (p.8|)), 
and the singularity is at t = — oo. The classical singularity can be avoided only if we set po = 0. But then, due to 
Equation ( |2.6| ), we would also have pi = 0, which corresponds to the trivial case of flat space-time. Therefore, the 
unique classical solution which is non-singular is the trivial flat space-time solution. 

III. QUANTIZATION OF THE CLASSICAL MODEL 

The Dirac quantization procedure yields the Wheeler-DeWitt equation through the imposition of the condition 

7i:* = 0, (3.1) 



on the quantum states, with Ti. defined as in Equation (2.6) (we are assuming the covariant factor ordering) using the 
substitutions 



Equation (3.1) reads 



^fl^*^ " (3 - 3) 
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A. The empty model 



In a first and simple example, we will freeze two degrees of freedom: the matter degree of freedom y 3 = <ft = 0, and 
one of the anisotropic degrees of freedom, y 2 = (3- — 0. We obtain a two-dimensional Klein-Gordon equation whose 
general solution is 

(3.4) 



(3.5) 



¥(#,, (3+) = J {F(k) exp[ik((3 + 13+)] + G(k) exp[ifc(/3 - (3+)}}dk 

with F(k) and G{k) arbitrary functions of k. Let us take gaussian superpositions 

(k-d) 2 ' 



F{k) = G(-k) = exp 



where a and d are constants. In this case the wave function reads 



((3 Q + f3+) 2 a 2 



\&1 = fT\/ff^ ex P 

For the case where 



the wave function turns out to be 



exp[id(/3 + (3+)] + exp 



(/3 -/3+)V 



exp[-irf(/3 - 13+)] 



F(k) = G(k) = exp 



(k + df 



* 2 = er%/7r jexp 



((3 + f3+) 2 a 2 



exp[-id(/3 + (3+)] + exp 



(/?o - (3+) 2 o 2 



exp[-id(/3 - /3+)] 



(3.6) 



(3.7) 



(3.8) 



In the next section we will calculate the bohmian trajectories corresponding to the solutions (3.6,3.8) and see if 
there are isotropic phases in these quantum models. 



B. The scalar field model 



In the g ener al case of the four-dimensional minisuperspace model, we will investigate spherical-wave solutions of 
Equation (3.3). They read 



<jr, = - 



f(y° + y) + g(y°-y) 



where y = yEj=iM 2 - 

One particular example is the gaussian superposition of plane wave solutions of Equation (3.3) 

= I {F{k) exp[z(|%° + k.y)] + G(fc) exp[z(|%° - k.y)]}d 3 k , 



where k = {ki,k 2 ,k 3 ), y = (y 1 ^ 2 ^ 3 ), |fe| = Jj2 l= i( k i) 2 > with F ( k ) and G ( k ) S iven b Y 



F(k) = G{k) = exp 



(\k\-d) 2 



After performing the integration in Equation ( 3.10| ) using spherical coordinates we obtain (2 



r 3/2 



*4(y°,y) = [2dtr + i{y° - y)<i 3 ] exp 



exp[id(y° - y)] 



i •['!-■ ,i,r y)- 2 



-[2da + i{y" + y)a 3 ] exp 



(y" + y) 2 a 2 



exp[id(y° + y)] 



i - <l»l - + t(y° + y)^ 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



where $(x) = (2/y/n) J Q exp(—t 2 )dt is the probability integral. The wave function ^4 is a spherical solution with the 



form of E quatio n (3.9) with g = — f. In order to simplify $4, we will take the limit a 1 >> d and (y° ± y)a >> 1 in 



Equation ( p.l2[ ) yielding 



IV. THE QUANTUM BOHMIAN TRAJECTORIES AND RESULTS 

In this section, we will apply the rules of the Bohm-de Broglie interpretation to the wave functions we have obtained 
in the previous section. We first summarize these rules for the case of homogeneous minisuperspace models. In this 
case, the general minisuperspace Wheeler-DeWitt equation is 

H[p a (t),q a (m(q) = 0, (4.1) 

where p a (t) and q a (t) represent the homogeneous d egre es of freedom coming from the gravitational and matter fields. 



Writing $ = Rcxp(iS/h), and substituting it into (4.1), we obtain the equation 



1 f)Q f)Q 

where 

The quantities / a/3 (g M ) and U^) are the minisuperspace particularizations of the DeWitt metric [E2|, and of the 
scalar curvature density of the space- like hypersurfaces together with matter potential energies, respectively. The 
causal interpretation applied to quantum cosmology states that the trajectories q a (t) are real, independently of any 
observations. Equation (|4.2|) is the Hamilton- Jacobi equation for them, which is the classical one amended with the 
quantum potential term ( }f .3| ), responsible for the quantum effects. This suggests to define 

P " = W ' (44) 
where the momenta are related to the velocities in the usual way 

1 dq 13 

Pa = Upjj-Qf > ( 4 - 5 ) 

and f a /3 is the inverse of f a ^. To obtain the quantum trajectories, we have to solve the following system of first order 
differential equations, called the guidance relations: 



Equations (4.6) are invariant under time reparametrization. Hence, even at the quantum level, different choices of 
N(t) yield the same space-time geometry for a given non-classical solution q a (t). There is no problem of time in the 
causal interpretation of minisuperspace quantum cosmology^. 

For the minisuperspace we are investigating, the guidance relations in the gauge N = 12 exp(3yo) are (see Equations 



where S is the phase of the wave function. 



1 This is not the case, however, for the full superspace (see Reference Jlq]) 
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A. The empty model 



Taking the phase of the solution ( |3.6| ), and inserting it into Equations (4.7), remembering that y 2 — (3- 
have been suppressed, we get the following system of planar equations: 



and y 



3 _ 



A) 



/3+cr 2 sin(2d/3 ) + 2dsinh(/3 /3+cr 2 
2[cos(2d/3 ) + cosh(/3 /3+cr 2 )] 



(4. 



2dcos(2d/3 ) + 2dcosh(/3 /3 + er 2 ) - /3 a 2 sin(2d/3 ) 
2[cos(2d/3 ) + cosh(/3 /3+cr 2 )] 



(4.9) 



The line (3q — divides configuration space in two symmetric regions. The line (3+ = contains all singular points of 
this system, which are nodes and centers. The nodes appear when the denominator of the above equations, which is 
proportional to the absolute value of the wave function, is zero. No trajectory can pass through them. They happen 
when (3 + = and cos(2g?/3 ) = — 1, or (3q = (2n+ l)ir/2d, n an integer, with separation n/d. The center points appear 
when the numerators are zero. They are given by /3+ — and (3q — 2d [cotan(d/?o )]/£>" . They are intercalated with 
the node points. As | (3q | — > oo these points tend to mr /d, and their separations cannot exceed n/d. As one can see 
from Equations (4.8) and (4.9), the classical solutions (vM are reco vered when | /3q |— > oo and | /3+ |— > oo, the /3's 



become proportional to t. The quantum potential given in Equation (4.3), which now reads 



Q(/3 , /?+) = -_ 



1 / d 2 R d 2 R 



df3 2 



dP 2 



[fl+o 2 sin(2d/? ) + 2dsinh(/3 /3 + cr 2 )] 2 - [/V 2 sin(2d/3 ) - 2dcos(2d/3 ) - 2dcosh(/3 /3+cr 2 )] 2 

8[cos(2d/3 ) + cosh(/5 0/ 3+cr 2 )] 2 



(4.10) 



becomes negligible in these limits. 
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FIG. 1. Field plot of the system of planar equations (4.5 -4^) for a = d = 1, which uses the Bohm-de Broglie interpretation 
with the wave function vPi, Equation (3.6). Each arrow of the vector field is shaded according to its true length, black 



representing short vectors and white, long ones. The four shades of gray show the regions where the vector field is pointing 
to northeast, northwest, southeast or southwest. The black curves are the nullcline curves that separate these regions. The 
trajectories are the white curves with direction arrows. 
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A field plot of this planar system is shown in Figure 1 for a = d = 1. We can see plenty of different possibilities, 
depending on the initial conditions. Near the center points we can have oscillating universes without singularities and 
with periodic isotropic phases, with amplitude of oscillation of order 1. For negative values of /3q, the universe arises 
classically from a singularity but quantum effects become important forcing it to recollapse to another singularity, 
recovering classical behaviour near it. Isotropic phases may happen near their maximum size. For positive values 
of (3q, the universe contracts classically but when (3$ is small enough quantum effects become important creating an 
inflationary phase which avoids the singularity. The universe contracts to a minimum size and after reaching this 
point it expands forever, recovering the classical limit when (3$ becomes sufficiently large. In this case isotropic phases 
may happen when the universe is near its minimum size. In both cases we see that isotropic phases happen only when 
quantum effects are important. We can see that, for /3q negative, we have classical limit for small scale factor while 
for p positive we have classical limit for big scale factor. In these models with d = 1, the isotropic phases last very 
shortly. Only for | d \ « 1 can the isotropic phases be arbitrarily large because, as said above, the separation of the 
singular points goes like ir/d. 

For the wave function the analysis goes in the same way but we have to interchange (3q with (3 + in Figure 
1. In this case, we have also periodic solutions but the others are anisotropic universes arising classically from a 
singularity, experiencing quantum effects in the middle of their expansion when they become approximately isotropic, 
and recovering their classical behaviour for large values of (3q. Depending on the initial conditions, the isotropic phase 
can be arbitrarily large. There are no further possibilities. 

B. The scalar field model 

Let us study the spherical wave solutions ( 3^ ) of Equation ( ) . The guidance relations ( f4.7| ) are 



' — (^) (4.11. 



p i = a i S = $S(?§±)=v*, (4.12) 



*3 

where S is the phase of the wave function. In terms of / and g the above equations read 



V 1 f(y° + y) + 9 (y°-y) '» ( j 



. i = y^(fy + y)- 9 '(y°-y) u (414) 



y V f(y° + y) + g(v° -v) 

where the prime means derivative with respect to the argument of the functions / and g, and is the imaginary 
part of the complex number z. 



From Equations (4.14) we obtain that 

dy J y J 

which implies that y l (t) = Ay°(t), with no sum in j, where the are real constants, c* = 1/cf, and c\ — c\ — 
C3 = 1. Hence, apart some positive multiplicative c onsta nt, kn owing about one of the y % means knowing about all 
y % . Consequently, we can reduce the four equations ( 4 . 1 3| ) and ( 4.14 ) to a planar system by writing y = C|y 3 |, with 



C > 1, and working only with y° and y 3 , say. The planar system now reads 

f'(y + C\y 3 \) + g'(y°-C\y 3 \ 
f(yO + C \y3\)+g(y°-C\y3\) 



.3 _ sign(2, 3 ) Jf{y a + C|y 3 |) - g'(y° - C\y 3 \) 



V C { f (yO + C \y3\)+g(y°-C\y3\) ' 
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Note that if / = g, y 3 stabilizes at y 3 = because y 3 as well as all other time derivatives of y 3 are zero at this line. 
As y l (t) = c*-y J (i), all y l {t) become zero, and the cosmological model isotropizes forever once y 3 reaches this line. 
Of course one can find solutions where y 3 never reaches this line, but in this case there must be some region where 
y = 0, which implies y 1 = 0, and this is an isotropic region. Consequently, quantum anisotropic cosmological models 
with / = g always have an isotropic phase, which can become permanent in m any cases. 

As a concrete example, l et u s take the gaussian $4 given in Equation ( |3.12| ). It is a spherical wave solution of the 
Wheeler-DeWitt equation (3.3) with / = —g, and hence it does not necessarily have isotropic phases as described 
above for the case / = g. Figure 2 shows a field plot of this system for d/a 2 = — 10~ 4 and C = 2. In the region 
\(f>\ » | /3o | we have periods of isotropic expansion because <f> — const, implies (3± — const.. Depending on the value 
of |0|, this isotropic phase can be arbitrarily long. These trajectories are periodic and without singularities. Due 
to quantum effects, no trajectory crosses the point /3q = <fi = 0. Figure 3 shows the bohmian trajectories for the 
wave function given in Equation (3.12), now with d/a 2 = 10~ 4 and C = 2. They have some qualitative behaviours 
different from the precedent case. The isotropic phases, now in contraction, are no longer periodic. They are parts of 
universes that contract anisotropically from infinity, experience a period of isotropy, and contract anisotropically to a 
singularity. Figure 3 shows that no trajectories crosses the point (3q = <fi = 0. 

In order to get some analytical insight over Figures 2 and 3, we present the planar system obtained from the 
guidance relations corresponding to the wave function ( |3.12 ) in the approximation ( 3. IS ) : 

m3f3l + 6C 2 /3U 2 -C7 4 4 ) , x 

a = _r . j..,. , (4.i8) 



p 2 M-c 2 ^ 2 



C 2 (j) 2 



16d o3 



aiiai-c 2 ^ 



2\2 



id? 



(3/3 2 + C 2 ^) 2 



(4.19) 



where we have reset y° = 0q and y 3 — <ft. This approximation is not reliable in the lines (3q = ±C<fi. As one can 
see immediately from these equations, (3± = const, whenever |0| >> |/3q|, and the sign of d defines the trajectories 
direction. 
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FIG. 2. Field plot of the system of planar equations ( 1.16| - 4~17] ) for a = 100, d = — 1 and C = 2, which uses the Bohm-de 
Broglie interpretation with the wave function ^4 given in Equation (3.12). Each arrow of the vector field is shaded according 
to its true length, black representing short vectors and white, long ones. The dark shade of gray shows the regions where the 
derivative of the vector field points clockwise (the light shade of gray means the opposite), and this shading allows to see the 
regions with arbitrarily long periods of isotropic expansion corresponding to periodic universes without singularities. There are 
also black or white curves with direction arrows corresponding to universes contracting from infinity to a singularity. Due to 
quantum effects, the trajectories do not cross the origin. 
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FIG. 3. Field plot of the system of planar equations (4.16 ^17| ) for a — 100, d — 1 and C — 2, which uses the Bohm-de 
Broglie interpretation with the wave function ^4, Equation (3.12). Each arrow of the vector field is shaded according to its true 
length, black representing short vectors and white, long ones. The dark shade of gray shows the regions where the derivative 
of the vector field points clockwise (the light shade of gray means the opposite) and this shading allows to see the regions with 
arbitrarily long periods of isotropic contraction. The trajectories are the black or white curves with direction arrows; they do 
not cross the origin due to quantum effects. 



V. CONCLUSION 



Adopting the Bohm-de Broglie interpretation of quantum cosmology, we have shown that quantum effects can 
generate an efficient alternative mechanism for the isotropization of cosmological models. Anisotropic classical models 
which never isotropize may present arbitrarily large isotropic phases during the course of their evolutions if quantum 
effects are taken into account, without needing to introduce any classical inflationary phase. The models studied were 
Bianchi I models in empty space or filled with a free massless scalar field. 

There are questions and developments which should be undertaken within this approach: the dependence of the 
above results on boundary conditions of the Wheeler-DeWitt equation, their generalizations to other Bianchi models, 
and the conditions for which quantum effects can also induce homogeneity on cosmological models which are clas- 
sically inhomogeneous. The last issue is by far the most interesting but also the most complicated one because the 
implementation of the Bohm-de Broglie interpretation for inhomogeneous cosmological models is much more subtle 
and involved Jig], These questions will be the subject of our future investigations. 
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